پژوهش‌های جغرافیای طبیعی. 8593 ۰۵۲ شسمارةٌ F‏ زمستان ۱۳۹۹ 
ص ۵۸۷- ۵۶۷ 


مقايسة عملکرد مدل‌های سری زمانی 54181314 و Holt-Winters‏ با روش‌های ھوش 
مصنوعی در بيش بينى طوفان‌های گردوغبار (مطالعة موردی: استان سیستان و بلوچستان) 
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جکیده 

هدف از این پژوهش مقايسة عملکرد مدل‌های سری زمانی Holt-Winters  SARIMA‏ با روش‌های هوشس 
مصنوعی شامل شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ و سیستم استنباط عصبی - فازی تطبيقى 
(ANFIS)‏ به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ در فصل آتی است. بدین منظورء 
از داده‌های ساعتی گردوغبار و کدهای سازمان جهانی هواشناسی در ينج ایستگاه سینوبتیک استان سیستان و 
بلوجستان با طول 5555 آماری ۲۵ساله (۲۰۱۴-۱۹۹۰) استفاده شد. نتایج نشان داد روش ۸۸۳15 نسبت به ساير 
روش‌هاء بهترین عملکرد را داشت و معیارهای ارزیابی MAE RMSE R‏ و NS‏ آن به‌ترتیب از ۰/۷۲ ۰/۵۷ 
۲۳ ۱ تا ۰/۹۵ ۰۰۱۵۱ 9۰۰/۴۰ ۰/۹۶ متغیر بود. همچنین. با افزایش شاخص متوسط FDSD‏ در 
ایستگاه‌ها (از ۱/۰۶ تا ۷/۱۱ دقت بيش بينى همة روش‌ها افزایش داست. بر همین اساس» در سری زمانی 
۸ ضريب همبستگی بين مقادیر مشاهداتی و پیش‌بینی‌شده شاخص 7959 از ۰/۶۴ به ۰/۷۹ افزایش 
یافت. برای روش‌های RBF dolt-Winters‏ و ANFIS‏ مقدار نیز ضریب همبستگی ias‏ تیب از ۰/۷۰ تا e [AV‏ 


۵۹ تا ۰/۹۲ و ۰/۷۲ تا ۰/۹۵ متغیر بود. درمجموع. با مقايسة روش‌های مورد استفاده. روش هوش مصنوعی 
5 بهترین و مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ بدترین عملکرد را داستند. 


وازگا نکلیدی: تحلیل سری زمانى» شبکه‌های عصبی, گردوغبار» ۸۷و PACF‏ نواحی بحران یگردوغباری. 


مه 


مشد مه 

يديدة گردوغبار یکی از مخاطره‌های طبیعی است و سالانه خسارت‌های زیادی در جنوب غرب و جنوب شرق کشور ایجاد 
می کند. گردوغبار یکی از منابع آلوده‌کنندة هوای پاک و عاری از آلودگی است و علت افزايش طوفان‌های گردوغبار در 
ایران به دلیل موقعيت جغرافیایی کشور با توجه به قرارگیری در کمربند خشک و نیمه‌خشک جهان است و از سوی 2592 
نزدیک‌بودن با کشورهایی همچون عربستان» عراق» و سوریه که دارای سطح وسیعی بیابان به دلیل عامل‌های مختلفی 
است. مانند کمبود پوشش گیاهی. بارش «Sul‏ و دمای بالا (عبدالشاه‌نژاد و همکاران» ۱۳۹۹: (VY‏ همچنین» وقوع 
خشک‌سالی در یک منطقه با تأمين آب در أن منطقه رابطة تنگاتنگی دارد و هر جه اقليم منطقه خشک‌تر باشد» احتمال 
وقوع خشکی و خشک‌سالی بیشتر شده و رخداد يديدة گردوغبار محتمل‌تر است (شاکر سوره و اسدی» (Ae AA‏ 


' نويسندة مسئول تلفن: ۰۹۱۲۲۱۵۹۷۳۸ Email: AliaghatQut.ac.ir‏ 


به طور كلىء منابع عمدة گردوغبار عبارت‌اند از: سرزمین‌های كمارتفاع و کمبارانی که از ميانكين بارش سالیانه کمتر 
از ۲۵۰ میلی‌متر برخوردارند. کشور ایران نیز از این قاعده مستثنى نیست (گودی» ۲۰۱۳: ۱۰۱). ايران» به ee‏ فرارگرفتن 
در کمربند خشک و نیمه‌خشک» در معرض سیستم‌های متعدد گردوغبار محلی و فرامنطقه‌ای است. همچنین, نيمة غربی 
کشور به دليل عرصه‌های بیابانی و بادهای منطقه‌ای» هرساله آسیب‌های زیادی از اين مخاطره را تجربه می کند (حسین 
حمزه و همکاران» ۱۳۹۵: [AY‏ 

در بين متغیرهای اقلیمی» بارش و دما از اهمیت ب‌سزایی برخوردارند و شبیه‌سازی دقیق بارش و دما و ارتباط أن با 
گردوغبار در مطالعات هواشناسی و کشاورزی از اهمیت زیادی برخوردار است. با وجود شبیه‌سازی بارش و دما به دليل 
وجود همبستگی درونی بين آن‌ها مشکل است (تانارته و همکاران» ۲۰۱۲). نتایج علمی مشابه توسط هرویجر و همکاران 
(۲۰۰۷) و هاننبرگر و نیکول (۲۰۱۴) نشان داد با افزايش دما تعداد وقوع طوفان‌های گردوغبار نیز افزايش مىيابد. 
همچنین, نتایج اين تحقیق در ارتباط با همبستگی پارامترهای اقلیمی با طوفان‌های گردوغبار همسوست با نتایج مطالعات 
sj‏ و همکاران (۲۰۰۶) و آمگالان و همکاران (۲۰۱۷)» که افزایش دماء افزايش سرعت ob‏ و کاهش بارندگی را از 
مهم‌ترین دلایل افزايش رخداد پديدة گردوغبار معرفی کرده‌اند. آیگنگ و همکاران (۲۰۰۹) با بررسی ارتباط گرمایش 
جهانی و تغییرات طوفان‌های گردوغباری جينء به اين نتيجه رسیدند که گرمایش جهان سبب کاهش شيب عرضی دما و 
در نتيجه کاهش شدت باد می‌شود. 

جمالی‌زاده تاج‌آبادی و همکاران (۱۳۸۶) توانایی دو روش شبکه‌های عصبی مصنوعی (ANN)‏ و ماشین‌های بردار 
پشتیبان (SVM)‏ در پیش‌بینی طوفان‌های گردوغبار شهر زابل را با طول دورة آماری ۲۶ساله (۲۰۰۵-۱۹۸۰) بررسی 
کردند. از exl‏ كه استفاده از نوع 7-5۷ و تابع کرنل Ab‏ شعاعی (RBF)‏ بهترین جواب‌ها را داده است» توصیه شده در 
مطالعات دیگر نیز قابلیت‌های اين روش‌ها بررسی شود. 

سبحانی و همکاران (۱۳۹۴) به بررسی گردوغبار و ارزیابی امکان پیش‌بینی أن بر اساس روش‌های آماری و مدل 
۳ در زابل با آمار چهارده‌ساله پرداختند. نتايج پیش‌بینی گردوغبار با مدل ANFIS‏ نشان‌دهندة cub‏ بالای أن در 
پیش‌بینی گردوغبار در اين ایستگاه است. علی‌یاری و همکاران (۲۰۰۸) نیز به پیش‌بینی کوتاه‌مدت آلودگی هوا به کمک 
شبکه‌های عصبی پرسپترون چندلایه J(MLP)‏ خط حافظه‌دار تأخير (DLM)‏ گاما و سیستم استنتاج فازی- عصبی 
تطبیقی (ANFIS)‏ پرداختند. نتيجه نشان داد روش پیشنهادی ترکیبی بر اساس بهینه‌سازی ازدحام ذرات (PSO)‏ و فیلتر 
کالمن برای آموزش ANFIS AS,‏ توانایی مناسبی در بهبود عملکرد پیش‌بینی دارد. 

سروستان و فلاح قالهرى (YA)‏ به پیش‌بینی يديدة گردوغبار در هشت ایستگاه منتخب استان خوزستان با استفاده 
از مدل‌های سری زمانی باکس جنکینز پرداختند. amus‏ نشان داد الگوهای مناسب ماهانه به‌ترتیب برای رامهرمز 
ARIMA (Y e ANAN)‏ آغاجاری (۰۱۰۱()۱۰۱۰۱ ARMA (Y‏ بهبهان ARIMA (ENFAN)‏ آبادان (Y V)‏ 
۸ دزفول ARIMA (MY) y)‏ امیدیه (۸۳(۸۰۸) «ARIMA‏ اهواز (۰۸۳(۸۸) ARIMA‏ و 
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مقايسة عملکرد مدلهاى سری زمانی Holt-Winters  SARIMA‏ با روش‌های هوس مصنوعی در پیش‌بینی e‏ ۵۶۹ 


مسجدسلیمان (۰۳۰۴()۱۰۱۰۱) ARIMA‏ می‌باشند که اين مدل‌ها از دقت بسيار خوبى برای پیش‌بینی گردوغبار برخوردار 
نبودند و همجنين بيش بينى تعداد روزهای طوفان گردوغبار برای سال‌های ۲۰۱۸ تا ۲۰۲۷ نشان داد که آغاجاری, «oll‏ 
و مسجدسلیمان بیشتر با پديدة گردوغبار مواجه‌اند. 

گودرزی و روزبهانی (۱۳۹۶) کارایی مدل‌های سری زمانی ARIMA‏ و هالت وینترز در پیش‌بینی دما و بارش ماهانة 
ایستگاه لتیان را بررسی کردند. gelo‏ پیش‌بینی توسط مدل‌های ARIMA‏ فصلی نشان از قابلیت بالای اين مدل‌ها در 
پیش‌بینی دمای ماهانه داشت و مدل ۲( ۰۰۱۲ (۰.۰۰۱) ARIMA‏ به‌عنوان مدل سری زمانی مناسب برای داده‌های دما 
تشخیص داده شد. اگرچه مدل هالت- وینترز نيز از دقت بالايى در پیش‌بینی دما برخوردار بوده نسبت به مدل ARIMA‏ 
خطای بیشتری داشت. همجنین» نتایج نشان داد مدل‌های ARIMA‏ فصلی و مدل هالت-وینترز قابلیت بالایی در 
پیش‌بینی مقادیر بالای بارش ماهانه در ایستگاه لتیان ندارند. اين مدل‌ها عمدتاً مقادیر بارش‌های پایین و متوسط را بهتر 
برآورد می کنند و تخمینگر مناسبی در بارش‌های بالا نیستند. 

سبحانی و همکاران (Y)‏ به بررسی مدل‌سازی و پیش‌بینی گردوغبار در غرب ايران پرداختند. بدین منظور از داده‌های 
گردوغبا do»‏ و رطوبت YA‏ ایستگاه مناطق درگیر شدید با گردوغبار در ايران در 85b‏ زمانی ٩۲ساله‏ (۲۰۱۸-۱۹۹۰) استفاده شد. 
بر اساس نتایج پژوهش. در بهترین شرایط مقدار RMSE‏ مدل ANFIS‏ برابر با ۱۱,۶۷ و مدل RBF‏ برابر با ۲,۱۹ است. بنابراین» 
قدرت RBF cà»‏ در پیش‌بینی گردوغبار در سال‌های شبیه‌سازی‌شده بیشتر است. بر اساس نتایج خروجی مدل RBF‏ در 
پیش‌بینی گردوغبار برای سال‌های آتی ایستگاه‌های مورد مطالعه. در هر دو مقیاس میانگین و حداکثر فراوانی كردوغبار, 
ایستگاه‌های غربی و جنوب غربی منطقة مورد مطالعه بیشتر در معرض گردوغبار در سال‌های أينده قرار گرفتند. 

انتظاری و سروستان (۱۳۹۶) به بررسی گردوغبار و پیش‌بینی أن در استان خوزستان برای باز زمانی ۲۰۲۳-۲۰۱۹ با 
a an IST‏ ی اساسا دای ار باه شام مات ارس 
خوزستان در طول دورة آماری ۲۰۱۰-۱۹۹۰ استفاده شد. نتیجة یافته‌های مناسب‌ترین مدل‌های سری زمانی گردوغبار 
برای شهرهای مسجدسلیمان و بهبهان مدل پیش‌بینی هالت-وينترزء شهر دزفول به‌تنهایی از مدل ARIMA (Y ۰۰ Y)‏ و 
شهرهای اهواز رامهرمز آغاجاری, آبادان» و اميديه دارای مدل نمو هموار ساده به‌دست آمد. 

تغییر اقلیم و بررسی اثرهای أن بر منابع آب از موضوعات مهمی به‌شمار می‌رود که پژوهشگران زیادی هدف تحقیق 
خود را در اين زمينه تعريف کرده‌اند. بنابراین» آهمیت پیش‌بینی وقوع پدیده‌های جوی نظیر كردوغبارء که بسیار حاثر 
همیت است. بیش از پیش حس می‌شود. تأثیر يديدة گردوغبار در ایران أن قدر وسیع است که بیش از نیمی از 
استان‌های کشور را به نحوی با مسائل و محدودیت‌های اين يديدة طبیعی درگیر کرده است که» علاوه بر اثرهای 
زیست‌محیطی» موجب اختلال در اجرای طرح‌های توسعة پایدار ملی شده و تاکنون پیامدهای منفی زیادی به دنبال 
داشته و خواهد داشت. از طرفی» بیشتر مطالعات داخلی در اين زمينه مربوط به روند پدیده‌های گردوغبار در محدوده‌های 
کوچک. مطالمات سینوپتیکی, و ماهوارای آن است. culla‏ با توجه به اينکه این پدیده اثرهای تامطلوب و پیامدهای 
منفی در زمینه‌های اجتماعی» اقتصادی» و سلامت مردم داشته است. مطالعه و پیش‌بینی و ارتباطسنحی أن با نوسانات 
اقلیمی ضروری است. در اين مطالعه به مقايسة عملکرد مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ با روش‌های 


هوش مصنوعی RBF‏ و ۸۱۷۳5 به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ پرداخته شده 


]. Frequency of Dust Stromy Days 


است و بنابراین اين تحقیق می‌تواند راهنما و الگوی جدیدی به‌منظور پیش‌بینی يديدة طوفان گردوغبار در کشور باشد؛ 
زيرا تا کنون چنین مطالعه‌ای با اهداف فوق‌الذ کر انجام نشده است. 


مواد و روس‌ها 
منطقة مورد مطالعه و داده‌های مورد استفاده 


در اين پژوهش به بررسی مقايسة عملکرد روش‌های آماری کلاسیک SARIMA)‏ و (Holt-Winters‏ و مدل‌های هوش 
مصنوعی RBF)‏ و (ANFIS‏ به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ در ينج ایستگاه 
سینوپتیک استان سیستان و بلوچستان (زابل زاهدان, ايرانشهرء خاش, و سراوان» که بیشترین فراوانی روزهای همراه با 
طوفان‌های گردوغبار در کشور را شامل می‌شوند (عراقی‌نژاد و همکاران؛ ۱۳۹۷: ۲۱ با جامعة آماری ۲۵ساله (۱۹۹۰- 
۴ در مقیاس فصلی يرداخته شد. بدین منظور از داده‌های ساعتی قدرت دید افقی و کدهای سازمان جهانی 
هواشناسی (WMO)‏ استفاده شد. مشاهدات پدیده‌های هواشناسی به فاصلة سه ساعت یک بار و درمجموع هشت بار در 
شبانه‌روز ثبت می‌شود. در اين مشاهدات» پدیده‌های بصری أب و هوا بر اساس دستورالعمل سازمان جهانی هواشناسی در 
اک که TEENS E‏ کی ET R1 SO E‏ کار کر انمق قفا 
مختلف هواشناسی از ۱۱ كد استفاده می‌شود (اولینگ‌سای و همکاران» ۲۰۱۴). کدهای مربوط به طوفان‌های گردوغبار در 
جدول اارائه شده است. با توجه به اينكه داده‌های گردوغبار به‌صورت تفکیک‌نشده از سایر پدیده‌های هواشناسی از 
سازمان هواشناسی در اختیار قرار می‌گیرد» در وهلة اول کدهای مربوط به گردوغبار از ساير پدیده‌های اقلیمی در ساعات 
مختلف همدیدی جدا شد (گودی و میدلتون. ۲۰۰۶). بر اساس تعریف سازمان جهانی هواشناسی روزهای همراه با 
طوفان‌های گردوغبار به روزی گفته می‌شود که حداقل در یکی از هشت سینوپ (گزارش‌های سه‌ساعتة ديدهبانى) یکی از 
کدهای مربوط به گردوغبار (۰۶» ۰۷ ۰۸ ۰۹ ۲۰ تا YO‏ و (AA‏ در بخش هوای حاضر گزارش شده باشد. البته, به شرط 
اينكه داده‌های قدرت دید افقی متناظر با آن كد گردوغبار کمتر از ۱۰۰۰ متر به ثبت رسیده باشد (گودی و میدلتون» 
EIE nua‏ شیرتا ماک رفن ایا و TEM‏ ي 
هواشناسی گردوغبار استفاده شده است. 


جدول ۱. کدهای سازمان جهانی هواشناسی مرتبط با فرسایش بادی و پدیده‌های گردوغبار (اولینگ‌سای و همکاران. ۲۰۱۴) 


مه ناشی از گردوغبار 

id‏ غبار يا شن برخاسته از زمين 

۸ طوفان كردوغبار 

v‏ طوفان گردوغبار اتفاق‌افتاده در گذشته (اتفاق‌افتاده در یک ساعت قبل از مشاهده يا در ایستگاه) 

۳۰ گردوغبار خفیف يا متوسط شن و ماسه همراه b‏ کاهش میدان دید کمتر از ٠٠٠١‏ متر اما بیشتر از ۲۰۰ متر 

Y‏ گردوغبار پایدار يا متوسط شن و ماسه همراه با کاهش میدان دید کمتر از ٠٠٠١‏ متر اما بیشتر از ۲۰۰ متر 

۳۲ شروع با افزايش گردوغبار خفیف يا متوسط شن و ماسه همراه b‏ کاهش میدان دید کمتر از ٠٠٠١‏ متر اما بیشتر از ۲۰۰ متر 
۳۴ طوفان گردوغبار پایدار شدید همراه با کاهش دید کمتر از ۲۰۰ متر 

۵ شروع يا افزايش طوفان گردوغبار شدید همراه با کاهش دید کمتر از ۲۰۰ متر 

۸ رعدوبرق با گردوغبار و يا طوفان شن و ماسه 


1. World Meteorological Organization 
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يس از انتخاب ایستگاه‌ها و بررسی داده‌ها در بازة زمانی ۲۵ساله (۲۰۱۴-۱۹۹۰) تعداد روزهای همراه با طوفان 
گردوغبار (FDSD)‏ برای ينج ایستگاه هواشناسی مورد مطالعه در استان سیستان و بلوچستان با استفاده از داده‌های دید 
افقی و کدهای سازمان هواشناسی محاسبه و در جدول Y‏ نشان داده شده است. ضمن اینکه ایستگاه‌های هواشناسی. 
طول و عرض جغرافيايى» ارتفاع از سطح درياء متوسط شاخص FDSD‏ در مقیاس فصلی و تعداد روزهای گردوغباری 
به‌ترتیب نزولی قابل مشاهده است. ایستگاه زابل» که در بين ابستگاه‌های مورد مطالعه در مرزی‌ترین نقطة استان سیستان 
و بلوچستان واقع شده با ۷۱۱ روز در بازة زمانی ۲۵ساله به‌عنوان رکورددار بیشترین روزهای همراه با طوفان گردوغبار در 
اين استان ثبت شده است. بر اين اساس» ایستگاه‌های خاش 9 سراوان به‌ثر تیب b‏ ۱/۸۰ 9 ۱۰۶ روز کمترین فراوانی 
روزهاى همراه با طوفان گردوغبار را در بين ابستگاه‌های مورد مطالعه به خود اختصاص دادهاند. بر اساس اين توضیحات 
به صورت جغرافیایی می‌توان بیان کرد که با پیشروی از مناطق شرقی و شمالی استان (مناطق مرزی) به سمت جنوب و 
جنوب غرب استان از وقوع فراوانی روزهای همراه با طوفان گردوغبار کاسته شده است. ede‏ أن را می‌توان در وزش 
بادهای ۱۲۰روزه» كه در نتيجة تشکیل سیستم پرفشار در شمال شرق دریای خزر و شکل‌گیری یک سامانة کم‌فشار در 
خشک‌شدن تالاب‌های منطقه جست‌وجو کرد. نتایج این بخش b‏ مطالعات مپرابی و همکاران (WAY)‏ انتظاری و 
سروستان (YAT)‏ عراقی‌نژاد و همکاران (۱۳۹۷). و انصاری قوجقار و همکاران (۱۳۹۹) همخوانی دارد. شکل ۱ موقعیت 


جدول ۲. فراوانی روزهای همراه با طوفان‌های گردوغبار استان سیستان و بلوچستان (۲۰۱۴-۱۹۹۰) 
عرض جغرافیایی طول جغرافیایی میانکین فصلی FDSD FDSD‏ 


ce‏ (درجه) ibd au ME‏ روز روز 
y ۱۹۵ 2-7 yyy- EUR‏ ۱.۶ 
خاش Me vis WAY BAM YAY‏ 
Yvy PEOR‏ ۲ع ۵۹۱ YM YA‏ 
زاهدان ۳۹۲۳۸ ۶۰۵۳ .۱۳۷ YA‏ ۳۴ 
زابل YM ۴۳۸۹ FAY yyy-‏ ۷/۱ 
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آنالیز اولیه داده‌ها 

برای آنالیزه مدل‌سازی» و پیش‌بینی از نرم‌افزارهای R‏ و Matlab‏ استفاده شد. پیش از هرگونه اقدام در جهت مدل‌سازی يا 
پیش‌بینی سری زمانی» تجزیه و تحلیل توصیفی اولیه‌ای از داده‌ها امری ضروری است. در اين مطالعه آنالیز اولية داده‌ها بر 
روی مجموعه داده‌های اموزشی صورت گرفت که شامل شناسایی ویژگی‌های مهمی نظیر همبستگی.» الگوهای فصلی. 
تغییرات دوره‌ای تودرتو در الگوهای فصلىء روند» داده‌های پرت» و هرگونه نوسانات قابل توجه دیگر در اين سری‌هاست. 
هه ddl‏ ادوا nay bs‏ کو کور از مستا نون انش ات که usb‏ ان اسان اما چ مات فا گید 3 
قاری ترش di a‏ باق io odo E So odia o‏ 
پیش‌بینی فرض می‌کنند داده‌های سری زمانی ایستا می‌باشند؛ زیرا با این فرض پیچیدگی ریاضی مدل‌های برازش 
داده‌شده کاهش می‌یابد؛ حال آنکه اغلب در واقعیت چنین نیست. در اين مطالعه از آزمون ریاضی Dickey-Fuller‏ برای 
تشخیص ایستایی b‏ ناایستایی دادهها استفاده شد (كوكران» ۰۰۵۵ : (YY‏ آنگاه در مواردى كه سری زمانى ناايستا بود از 
تفاضل گیری " و تبدیلات توان برای حذف ناایستایی روند يا الگوهای فصلی استفاده شد. آنالیز اولي داده‌ها نیز از طریق 
تجزیه و تحلیل بصری نمودارهای زمانی شامل ACF‏ ۳۸۲۳ و محاسبات آمار توصیفی پایه انحام گرفت. برای 
شناسایی الگوهای موجود در سری‌های زمانی از تابع Decompose‏ در نرم‌افزار R‏ استفاده شد که اين تابع از روش 
میانگین متحرک برای تحزیة سری زمانی به اجزای فصلىء روند» و بخش نامنظم و تصادفی استفاده می کند. 


انتخاب و برازس مدل‌ها 

دارد. تحقیقات گسترده‌ای انحام شده است که نشان می‌دهند یک روش پیش‌بینی برتر برای همه موارد و در حالت جامع 
پیش‌بینی با صحت قابل قبول» بايد روش‌های پیش‌بینی مختلفی مقایسه شود و توانایی هر یک در پیش‌بینی مشاهدات 
نمونه‌های آینده ارزیابی شود. به دلیل وجود الگوهای فصلی در داده‌های ایستگاه‌های مورد مطالعه» مدل‌های پر کاربرد و 
دارای صحت بالاتر در مطالعات مشابه. كه قابلیت بيش بينى الگوهای فصلی را نیز داشتند. به‌عنوان مدل‌های پیش‌بینی 
سرى زمانى در اين مطالعه درنظر گرفته شده‌اند كه اين مدل‌های آماری عبارت‌اند از مدل‌های سری زمانی Holt-‏ 
Winters‏ و SARIMA"‏ (ابوكىسارفو و همکاران» ۲۰۱۵) و مدلهاى يادكيرى ماشين شامل دو روش RBF‏ و ANFIS‏ 


(سينك و ميشراء ۲۰۱۵). 


مدل SARIMA‏ 
مدل SARIMA‏ گسترش يافتة مدل ARIMA?‏ برای داده‌های فصلى است. منظور از سری‌های زمانى فصلی» سرىهايى 
اهوت كه 2525 ous‏ تم Delos jio‏ ات obs‏ وا قرام e. o‏ کی کم کا قمعا من سرض 
زمانی» که دارای الگوهای فصلی‌اند. نه‌تنها بررسی 3999 همبستگی بين مشاهدات فعلی و مشاهدات اخیر ضروری است 


]. Stationary 

2. Differencing 

3. Auto Correlation Function 

4. Partial Autocorrelation Function 

5. Seasonal Autoregressive Integrated Moving Average 
6. Average Autoregressive Integrated Moving 
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که بایستی وجود همبستگی ميان مشاهدات فعلی و مشاهدات قبلی در فصل‌های گذشته نیز بررسی شود که روشی برای 
شناسایی وجود يا فقدان الگوی فصلی در داده‌ها نیز هست. اين مدل فرض می‌کند که سری زمانی درنظر گرفته‌ شده خطی 
است و از یک توزیم آماری شناخته‌شدة خاص, مانند توزیع نرمال» تبعیت می‌کند. محبوبیت مدل SARIMA‏ عمدتاً به 
دليل انعطاف‌پذیری أن در نمايش ساده انواع مختلفی از سری‌های زمانی برای فرایند ساخت مدل بهینه است. برای اجرای 
مدل SARIMA‏ از الگوریتم ارائه‌شده توسط هایندمن و خندکار استفاده شده است (هایندمن و خندکان ۲۰۰۸). 


مدل Holt-Winters‏ 
e 6 ۰ 5 ^ 3 n‏ * ۰ 
اين مدل یکی از محبوب‌ترین روش‌های هموارسازی نمایی است. اين مدل بر اساس الگوهای فصلی و روند موجود در 
داده‌ها ایحاد می‌شود. محبوبیت اين روش به دلیل مزایایی از قبیل نیاز محاسباتی كمترء سادگی استفاده» و دقت پیش‌بینی 
این مدل به‌خصوص برای سری‌های دارای الگوی فصلی است. اين روش دو گونه معادله Additive)‏ و (Multiplicative‏ 
دارد. از آنجاکه نوسانات سری‌های زمانی ابن مطالعه تسيا ثابت بردتت از معاد Additive‏ استفاده شد که برای سر 
زمانی lx: t- 1, ..., n!‏ به صورت Ča] B)‏ 
X^, ml t hb, + S (Y)‏ 
كه ظ افق پیش‌ببنی zu‏ برای مثال» در این مطالعه كه LU.‏ برای فصل آتی (سه ماه آینده) است» Y b yl» h‏ )039 
و ساير يارامترها عبارت‌اند از: 
Level: /, za(y, — s, ,,) ۲ (d-a) I, , b, 4)‏ 

Trend: b, -p*(I, — 1, (+ )1- B, , (Y) 

Seasonality: s, =y(y, — Ll, — b, 4) +A- 9s, 4 


که B* ax‏ و y‏ پارامترهای هموارسازی‌اند و مقدار أن به صورت 1 > © > 0 1 > *8 > ۵ و © -1 > ۲ > 0 بوده و 
همچنین m‏ طول 8555 فصلی است. 


شبكة عصبی مبتنی بر توابع پایة شعاعی (RBF)‏ 

شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ دارای پاية ریاضیاتی بسیار قوی بر مبنای فرضية منظم‌سازی برای حل 
مسائل مشکل است. به طور کلی اين شبکه از سه لايه شامل à Y‏ ورودی» مخفی» و خروجی تشکیل شده است (چن و 
همکاران. ۱۹۱: ۲۰۳). در اين شبکه تابع انتقال گوسی در لابة مخفی استفاده می‌شود و در ALY‏ خروجی تابع انتقال 
خطی است. RBF‏ يايه و اساس نوع خاصی از شبکه‌های عصبی به نام شبکه‌های عصبی آماری است. در شکل Y‏ یک 
شبکة RBF‏ نشان داده شده است. همان طور که در شکل مشحص iul‏ نورون RBF‏ یک تابع گوسی FT‏ ورودى اين 
تابح فاصلة آقلیدسی بين هر ورودی به نورون با بردار مشخ صشدة هم‌اندازة بردار ورودی است. اين تابح گوسی از رابطة 


f(X,b)-e' r-|X, - X,|*0.8326/h (V) 


]. Hyndman 
2. Khandakar 
3. Exponential Smoothing 


در اين رابطه X,‏ ورودی شبکه با خروجی نامعلوم» ,ل ورودی مشاهداتی در زمان يا مکان sib‏ / پارامتری است که 
پهنای تابع گوسی را کنترل می کند. خروجی اين تابع بين صفر تا یک متغیر است. محاسب۵ خروجی ,۲ بر اساس متغیر 


مستقل X,‏ به صورت زیر به‌دست می‌آید: 
Y. = LW * f(X „,b) + Bias (Y)‏ 
در این رابطه LW‏ و Bias‏ به‌ترتیب وزن ماتریس ارتباطات بين AY‏ مخفی و LY‏ خروجی و ماتریس بایاس AY‏ 


rac c op 


E f CX, $ D) = e Vx — x, || 0.8326 AF 
| 227 | 
X3, ل‎ RE 




















RBF سماتیک‎ .Y JS 5 


سیستم استنباط عصبی - فازی تطبیقی (ANFIS)‏ 
شیستم. استباط =a‏ فازی تطبیقی ba (ANEIS)‏ نكن ور سال ۱۹۹۴ معرفی شد ANEIS ۱۹۹۳ scil)‏ 


شبیه به یک شبكة عصبی چندلابه است؛ با اين تفاوت که علاوه بر الگوریتم‌های بادگیری شبکة عصبی. از منطق فازی 
نیز بهره می‌گیرد. یک مدل ANFIS‏ از ينج SY‏ تشکیل شده است؛ اين ينج لابه به‌ترتیب عبارت‌اند از: AY‏ ورود 
اطلاعات à Y‏ محاسبة وزن قوانین فازی. LY‏ نرمال‌سازی وزن‌های قوانین به‌دست‌آمده, Y‏ محاسبة قوانین, LY‏ 
جمع‌بندی و خروجی شبکه. در اين تحقیق تابح عضویت درنظرگرفته‌شده تابح عضویت دوزنقه‌ای و الگوریتم آموزش 
شبکه روش هیبریدی درنظر گرفته شد. شکل ۲ شماتیکی از ANFIS‏ را نشان می‌دهد. 





شکل Y‏ شماتیک سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ 
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مدل‌های پیش‌بینی 

به‌منظور پیش‌بینی کمی فراوانی روزهای همراه با طوفان‌های گردوغبار (FDSD)‏ چهار مدل مختلف استفاده شد. سپس» 
بر اساس معیارهای ارزیابی بهترین مدل در هر بيش بينى انتخاب شد. متغیر تعداد فصل (های) گذشته برای پیش‌بینی 
فصل آتی به‌کار گرفته شد که می‌تواند شامل یک دوه سهء و چهار فصل قیل از متفیر فراوانی روزهای همراه ب 
طوفان‌های گردوغبارباشد. در aa‏ برای پیش‌بینی متفیر فراوانی روزهای همراه با طوفان‌های گردوغبار در فصل بعده 
از سری زمانی ۳5۲ با یک گام تأخیر استفاده شد. بدین صورت که برای پیش‌بینی در زمان 1+ از مقدار آن در زمان t‏ 
استفاده شد. در مدل ۲ پیش‌بینی در فصل بعد بر اساس مقادیر فراوانی روزهای همراه با طوفان‌های گردوغبار تا دو فصل 
قبل استفاده شد و به طور مشابه در مدل ۳ و ۴ پیش‌بینی‌ها بر اساس مقادیر متفیر مورد نظر تا سه و چهار ab‏ قبل 
صورت پذیرفت. روابط ۸-۵ زیر در تکمیل توضیحات بیان‌شده ارائه می‌شود. شکل ۶ مدل‌های مورداستفاده در انتخاب 
فصل (های) گذشته برای پیش‌بینی در فصل آتی را نشان می‌دهد. روش‌های بيش بينى به کاررفته در اين پژوهش شامل 
RBF Holt-Winters ۸‏ و ANFIS‏ است. 





FDSD, ب‎ f (FDSD,,) ۱ مدل‎ (à) 
FDSD, „=f (FDSD,,,, FDSD +) مدل۲‎ )۶( 
FDSD, „=f (FDSD,,, FDSD, ,,, FDSD 4) ۲ مدل‎ (Y) 
FDSD, „=f (FDSD,,,,FDSD, 4, FDSD, ,رو‎ FDSD,, 4) ۴ مدل‎ (A) 
"o 
(Y) مدل‎ 
UNIS SARIMA 
Holt-Winters 
(f) مدل‎ RBF 
© ANFIS 

سكل ۴. ساختارهای مورد استفاده در پیش‌بینی گردوغبار 

آماده‌سازی داده‌ها 


يس از بررسی و کنترل کیفی آمار ایستگاه‌ها و رفع نواقص آماری» همگنی داده‌ها با استفاده از آزمون Run Test‏ ارزیابی 
و تصادفی‌بودن داده‌ها در سطح اطمینان ۹۵درصد پذیرفته شد. در ادامه» با توجه به وجود سری زمانی متغير فراوانی 
روزهای همراه با طوفان‌های گردوغبار اقدام به پیش‌بینی أن با استفاده از روش‌های .Holt-Winters .SARIMA‏ 
RBF‏ و ANFIS‏ به‌منظور استفاده از روش‌های هوش مصنوعی شد. از ۸۰ درصد داده‌ها برای مرحلة آموزش و ۲۰ 


در صد برای آزمون استفاده شد. پارامترهای مناسب رو شهاى مختلف از طریق ازمون و خطا تعیین شد. همچنین. 


بدمنظور جلوگیری از بیش‌برازش از صحت‌سنجی متقاطع استفاده شد. در روش صحت‌سنجی متقاطع در هر بار اجرای 
مدل برای پیش‌بینی با یک پارامتر مشخصء یک يا تعدادی از داده‌های مجموعه آموزش كنار گذاشته می‌شوند (به‌عنوان 
داده‌های صحت‌سنحی) و مدل بر اساس داده‌های باقی‌مانده آموزش می‌بیند. سپس داده‌های کنارگذاشته‌شده به مجموعة 
آموزش برگردانده شده و داده يا مجموعة دیگری از داده‌های آموزش به‌عنوان دادة صحت‌سنحی از اين مسير جدا می‌شود 
و بار دیگر شبکه آموزش داده می‌شود. اين فرایند تا زمانی که همه داده‌های آموزش در قسمت صحت‌سنحی استفاده 
شوند» ادامه می‌یابد. يس از أن بر اساس میانگین خطای مدل در تخمین داده‌های صحت‌سنحی, پارامتر (های) مناسب 
شبکه عصبی تعیین می‌شوند. برای پیش‌بینی از چهار مدل معرفی‌شده فوق الذکر استفاده شد. مدل مناسب بر اساس 
معیارهای ارزیابی برای هر روش در مقیاس فصلی استخراج شد. پس از «ol‏ تحلیل نتایج بر پایه بهترین مدل برای هر 
يك از روش‌ها صورت پذیرفت. 


معیارهای ارزیابی 


e. 1 ۱ eus ۰ ۰ ۰ ۰ * مه $ همه ^ ا‎ m 
جذر میانگین‎ (R) برای ارزيابى مدل‌ها و مقايسة نتایج روش‌های مختلف از چهار معیار ارزیابی یعنی ضريب همبستگی‎ 
استفاده شد. اين معیارها به‌صورت‎ (NS) و نش‌ساتکلیف‎ (MAE) میانگین قدر مطلق خطا"‎ (RMSE) مربعات خطا"‎ 


زیر تعريف می‌سوند: 





در اين رابطه‌هاء 01 polio‏ مشاهداتی در گام زمانی eli‏ 11 مقادیر پیش‌بینی‌شده در گام زمانی لام 0 ميانكين مقادیر 
مشاهداتی. Í‏ ميانكين مقادیر پیش‌بینی‌شده» و n‏ تعداد داده‌هاست. هر مدل که دارای R‏ و NS‏ بیشتر و RMSE‏ و 


شکل ۵ تغییرات شاخص FDSD‏ را در پنج ايستكاه سينويتيك منتحب استان سیستان 9 بلوچستان از منظر فراوانی 
روزهای همراه b‏ طوفان گردوغبار نشان می‌دهد. بر اساس اين شكل» شاخص FDSD‏ در ایستگاه‌های سراوان» خاش 


1. Correlation 

2. Root Mean Square Error 
3. Mean Absolute Error 

4. Nash-Sutcliffe 


مقايسة عملکرد مدل‌های سری زمانی Holt-Winters  SARIMA‏ با روش‌های هوس مصنوعی در پیش‌بینی AYY e‏ 


افزايش تعداد روزهای گردوغباری در ایستگاه زابل از پراکندگی نوسانات کاسته دة 9 بك شدت ]4 ن افزوده شده انش 9 
مقادير اوج گردوغبار به صورت متمرکز در كنار یکدیگر مشاهده می‌شوند که حکایت از وقوع طوفان‌های گردوغبار 
يىدريى در اين ایستگاه از سال ۲۰۰۰ به بعد دارد. با توجه به شکل ۵ وجود فصول متوالی عاری از طوفان‌های 
گردوغبار انقطاع و داشتن یک سری زمانی منفصل, و به تبع أن تشکیل طوفان‌هایی با فاصلة زمانی بیشتر نسبت به 
ایستگاه زابل در سال‌های ۳ تا ۲۰۰۱ در ایستگاه‌های سراوان» خاش ایرانشه 9 زاهدان مشهود sa‏ بر اساس 
۲۰۱۲-۲۰۸ را در ایستگاه‌های مورد مطالعه تشخیص داد. نتایج مشايهى از ویژگی‌های شاخص FDSD‏ در مقیاس زمانی 
فصلی توسط زینالی (۱۳۹۵) عراقی‌نژاد و انصاری قوجقار (۱۳۹۶) و عراقی‌نژاد و همکاران (۱۳۹۷) ارائه شده است. 
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گام زمانی (فصلی) 


سکل ۵. سری زمانی ساخص FDSD‏ در مقیاس زمانی فصلی برای ایستگاه‌های مورد مطالعه 


اولين قدم در پیش‌بینی سری‌های زمانی» شناسایی ویژگی‌های مهمی نظیر همبستگی, الگوهای فصلی. تغییرات 
دوره‌ای در الگوهای فصلی. روند و ابستایی يا غیرایستایی داده‌هاست. وجود همبستگی در ميان داده‌های سری زمانى» 
امکان اجرای مدل‌سازی و پیش‌بینی مقادیر بعدی سری زمانی را فراهم می‌کند. در اين مطالعه. وجود همبستگی با رسم 
نمودارهای ACF‏ و PACF‏ برای هر ایستگاه با فاصلة اطمینان ۹۵ درصد (خطوط نقطه‌چین) بررسی شد (شکل ۶). اگر 
معنی‌داری با صفر اختلاف دارد. همان‌طور که در نمودارهای ACF‏ و PACF‏ ایستگاه‌های مورد مطالعه مشاهده می‌شود. 
وقفه‌های زمانی معنی‌دار : نشان‌دهندة وجود همبستگی بين مقادیر زمانی است كه اجرای مدل‌سازی 9 پیش ‌بینی مقادیر 
أينده (فصل آتی) شاخص ۴۵52 را برای هر ينج ایستگاه مورد مطالعه ممکن می‌سازد. گام مهم در تجزیه و تحلیل 
سری‌های زمانی بررسی ایستایی يا غیرایستایی میانگین سری‌های زمانی است. با توجه به شکل ۶ نمودارهای ACF‏ و 
PACF‏ ایستگاه‌ها حکابت از ایستایی داده‌های همه سری‌های زمانی به غير از ابستگاه زابل دارد. 


پژوهش‌های جغرافیای طبیعی» 5,93 ۵۲ شمارة F‏ زمستان ۱۳۹۹ 





Autocorrelation Function for saravan 
(with 596 significance limits for the autocorrelations) 


Autocorrelation 








Autocorrelation Function for khash 
(with 596 significance limits for the autocorrelations) 
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Autocorrelation Function for iranshahr 
(with 596 significance limits for the autocorrelations) 
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Autocorrelation Function for ۵۱ 
(with 5% significance limits for the autocorrelations) 




















5 
5 
5 
4 - 
Lag 
Autocorrelation Function for zabol 
(with 596 significance limits for the autocorrelations) 
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Partial Autocorrelation Function for saravan 
(with 5% significance limits for the partial autocorrelations) 


Partial Autocorrelation 

















Partial Autocorrelation Function for khash 
(with 5% significance limits for the partial autocorrelations) 





Partial Autocorrelation 
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Partial Autocorrelation Function for iranshahr 
(with 5% significance limits for the partial autocorrelations) 





Partial Autocorrelation 
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Partial Autocorrelation Function for zahedan 
(with 5% significance limits for the partial autocorrelations) 


Partial Autocorrelation 
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Partial Autocorrelation Function for zabol 
(with 596 significance limits for the partial autocorrelations) 


Partial Autocorrelation 
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شکل JP‏ نمودار خودهمبستگی (ACF)‏ و خودهمبستگی جزئی (PACF)‏ شاخص FDSD‏ در ایستگاه‌های مورد مطالعه 


مقايسة عملكرد مدلهاى سری زمانی Holt-Winters  SARIMA‏ با روش‌های هوس مصنوعی در پیش‌بینی ... ۵۹ 


در داده‌های غيرفصلى و سالانه 00-75 با يك بار تفاضل‌گیری درجة اول داده‌ها ایستا می‌شوند؛ ولی» با توجه به 
فصلی‌بودن و مقیاس داده‌هاء با تفاضل‌گیری درجة دوم (D22)‏ داده‌های ایستگاه زابل ایستایی شدند که نمودار توابع 
خودهمبستگی (ACF)‏ و خودهمبستگی جزتی (PACF)‏ أن بعد از تفاضل‌گیری در شکل V‏ نشان داده شده است. با توجه 
به توابع خودهمبستگی و خودهمبستگی جزئی» محدودة تغییرات آتورگرسیون و ميانكين متحرک تعیین و با استفاده از 
معیارهای ارزیابی مناسب بهترین مدل سری زمانی برای هر ایستگاه استخراج شد. همچنین, از آزمون Dickey-Fuller‏ 
نيز برای تشخیص ایستایی يا ناایستایی سری‌های زمانی بر روی مجموعه داده‌های آموزشی استفاده شد (جدول Y‏ در 
آزمون 1016-۳1167 سطح معنی‌داری برای احتمال ۰,۰۵ P-Value‏ درنظر گرفته شد. فرض این آزمون نالیستایی 
سری زمانی را بررسی می کند. بر اساس اين آزمون» صرفاً سری زمانی ایستگاه‌های زابل ناایستاست که نتايج نمودارهای 
ACF‏ و PACF‏ ایستگاه‌های مورد مطالعه را تأييد می کند. 





Autocorrelation Function for zabol.new Partial Autocorrelation Function for zabol.new 
(with 5% significance limits for the autocorrelations) (with 596 significance limits for the partial autocorrelations) 
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V JS‏ نمودار خودهمیستکی (ACF)‏ 9 خودهمیستکی جزئنی (PACF)‏ بعد از تفاضل کیری فصلی مرتبه دوم (D-2)‏ در ایستگاه زابل 


جدول .Y‏ نتایج آزمون Dickey-Fuller‏ 


ابستكاه P-Value‏ 
سراوان JA‏ 
خاش asy‏ 
ایرانشهر ۱ 
زاهدان 7“ 
زابل yA‏ 


جدول ۴ polio‏ معیارهای ارزیابی مربوط به مرحلة آزمایش روش‌های مختلف به‌منظور پیش‌بینی فراوانی روزهای 
همراه با طوفان گردوغبار را نشان می‌دهد. نتایج پیش‌بینی شاخص 170510 حاکی از عملکرد مناسب روش‌های سیستم 
استنباط عصبی- فازی تطبیقی (ANFIS)‏ و شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ در همة ایستگاه‌های مورد 
مطالعه (زابل» «ual;‏ ایرانشهر» خاش و سراوان) است. بر اساس جدول ۴ در ایستگاه‌های زابل» زاهدان» و ایرانشهر 
روش‌های سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ مدل‌های 
سری زمانی dHolt-Winters‏ و SARIMA‏ به‌ترتیب رتبۂ اول تا چهارم بهترین روش‌های بيش بينى شاخص FDSD‏ را بر 
اساس معیارهای ارزیابی به خود اختصاص دادند. برخلاف ایستگاه‌های مذكورء در ایستگاه خاش بعد از روش ۸۱۳۲5 


وو 


استفاده از مدل سری زمانی Holt- Winters‏ برای بيش بينى شاخص FDSD‏ عملکرد بهتری نسبت به روش ۸8۴ از خود 
نشان داد. به طوری كه ضرایب همبستگی (R)‏ و نش‌ساتکلیف (NS)‏ بين polio‏ مشاهداتی فراوانی روزهای همراه با 
طوفان گردوغبار و مقادیر پیش‌بینی‌شدة أن به‌ترتیب از ۰۷۱ به ۰۷۴ و ۰۷۰ به ۰,۷۱ و همچنین معیارهای خطای 
RMSE‏ و MAE‏ به‌ترتیب از ۰۸۷۰ به ۰,۶۹ و ۰,۵۶ به ۰,۵۵ بهبود يافت. همچنین» در ایستگاه سراوان نیز بعد از روش 
5 روش RBF‏ و مدل سری زمانی Holt-Winters‏ عملکرد تقریباً یکسانی از خود نشان دادند؛ به طوری که 
معیارهای ارزیابی MAE RMSE R‏ و NS‏ در روش RBF‏ نسبت à‏ روش Holt-Winters‏ به‌ترتیب از ۰,۶۹ ego‏ 


وو 


«RÀ 9 ۵‏ بك uS -7Y 9 ۵۶ A e.‏ یافت. b‏ نوجه بك اينكه معيارهاى ارزيابى RMSE‏ و NS‏ نشان دهندة 
برتری روش RBF‏ نسبت به مدل Holt-Winters‏ و معیارهای ارزیابی R۸‏ و MAE‏ حكايت از برتری مدل Holt-Winters‏ 
ایستگاه داشت. علت بهبود نتایج مدل سری زمانی Holt-Winters‏ در ایستگاه‌های خاش و سراوان را می‌توان در کاهش 
وقوع فراوانی روزهای همراه با طوفان گردوغبار جست‌وجو کرد؛ به طوری که اين موضوع حاکی از قابلیت بالای مدل 
سری زمانی Holt-Winters‏ در بيش بينى مقادیر کم شاخص FDSD‏ است؛ (Ja‏ با توجه به آهمیت موضوع پیش‌بینی 
طوفان‌های گردوغبار در مناطق بحرانی کشور (پرتعداد» نمی‌توان چندان به توانایی‌های اين مدل در پیش‌بینی فراوانی 
روزهای همراه با طوفان گردوغبار اتکا کرد. همجنينء نتایج نشان داد که مدل سری زمانی SARIMA‏ نسبت به ساير 
روش‌ها قابلیت (e Yo‏ در پیش‌بینی شاخص FDSD‏ در هیچ یک از ایستگاه‌های مورد مطالعه ندارد. نتایج اين بخش با 


جدول ۴. معیارهای ارزیابی محاسبه‌شدة مدل‌های آماری کلاسیک سری زمانی و روش‌های هوش مصنوعی برای پیش‌بینی شاخص FDSD‏ 


حالت معیارهای ارزیابی سراوان <١‏ خاش ايرانشهر زاهدان زابل 
مدل بهینه Y 3 y Y‏ ۳ 
(Y) ۰۸۷۰ )۴( FA (Y) FF (Y) R‏ اعم (Y)‏ ۸۷۹ 
(f) ۰,۷۹ )۴( ۰,۶۹ (F) NY )۴( ۷۵ (f) RMSE SARIMA‏ ۰,۶۵ 
(Y) ۰۸۵۷ (Y) ۰۸۵۶ (Y) ۰۸۵۷ (Y) ۰,۶۰ (Y) MAE‏ ۰۸۵۵ 
(F) FA (Y) ۰,۶۲ (Y) NS‏ اعم (Y) ۰,۶۰ (Y)‏ ۰۸۱ 
مدل بهینه Y Y y Y Y‏ 
AY (Y) ۰,۶۲ (Y) ۸۷۵ (V) ۸۷۴ (Y) ۸۷۰ (Y) R‏ 
(Y) ۰۸۷۸ (V) ۰,۶۷ (v) ۰,۶۹ (Y) ۰۸۷۲ (v) RMSE Holt-Winters‏ ۰,۶۲ 
(V) ۰۸۵۲ (V) ۰۸۵۲ (V) ۰۸۵۵ (Y) ۰۸۵۶ (Y) MAE‏ ۰۸۵۱ 
(V) ۰,۶۲ (v) ۰۸۷۰ (v) ۰۸۷۱ (Y) ۰,۶۷ (N) NS‏ ۰۸۵ 
مدل بهینه Y y F Y Y‏ 
«AY (Y) ۰,۶۷ (Y) ۰۸۷۶ (Y) ۸۷۱ (V) ۰,۶۹ (V) R‏ 
(Y) ۰۷۳ (Y) ۰,۶۵ (Y) ۰,۷۰ (Y) ۰,۷۱ (Y) RMSE RBF‏ ۰۶۰ 
(Y) ۰,۵۰ (Y) ۰/۴۹ (Y) ۰۸۵۶ (V) ۰۸۵۹ (V) MAE‏ ۰,۴۸ 
(Y) NS‏ اعم (۳) ۰۸۷۰ (Y) ۰,۶۵ (Y) ۰۸۷۱ (Y)‏ ۰/۱ 
مدل بهینه Y Y Y Y Y‏ 
(Y) ۰۸۷۸ (Y) ۸۷۷ (V) ۰۸۷۶ (V) ۸۷۲ (V) R‏ ۰/۵ 
(Y) ۰۵۲ (Y) AY (Y) AY () AY (Y) RMSE ANFIS‏ ۰۸۵۱ 
(V) ۰/۴۱ (Y) IY (V) ۰/۴۱ (Y) ۸۴۲ (V) MAE‏ ۰۸۴۰ 
AF (V) ۸۷۵ (V) ۸۷۲ (V) ۸۷۲ (V) ۰۸۷۱ (Y) NS‏ 
متوسط روزهای همراه با طوفان گردوغبار در مقياس فصلی Y" ۳۴ Y ۱۸ MPO‏ 
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با توجه به اينكه در جدول Y‏ ترتیب ایستگاه‌ها به لحاظ شاخص FDSD‏ از ایستگاه سراوان تا زابل به صورت صعودی 
odg‏ و متوسط شاخص FDSD‏ در مقیاس فصلی از ۱,۰۶ تا ۷,۸۱۱ متغیر است» می‌توان دریافت که عملکرد هم روش‌های 
پیش‌بینی شاخص ۳5۲ با افزايش تعداد روزهای همراه با طوفان گردوغبار رابطة مستقیمی دارد. بدین معنی که با 
افزايش فراوانی روزهای همراه با طوفان گردوغباره مقادیر همه معیارهای ارزیابی مورد استفاده بهبود می‌یابند؛ به نحوی 
که در روش SARIMA‏ ضریب همبستگی بين polio‏ مشاهداتی و بيش بينىشدة شاخص FDSD‏ از ۰/۶۶ به ۰۷۹ 
افزايش می‌یابد. همچنین» برای روش‌های RBF Holt-Winters‏ و ANFIS‏ نيز به‌ترتیب از ۰۷۰ تا CUN‏ ۰,۶۹ تا ۰/۲ 
و ۰/۷۲ تا ۰/۹۵ در ایستگاه‌های مورد مطالعه متغیر است. نتيجة دیگری که از جدول Y‏ قابل استنتاج است عملکرد بهتر و 
مناسب همة روش‌های پیش‌بینی شاخص ۳۲5۲ بر اساس معیارهای ارزیابی در ایستگاه ایرانشهر نسبت به ایستگاه 
زاهدان است؛ اين در حالی است که فراوانی روزهای همراه با طوفان‌های گردوغبار در ایستگاه ایرانشهر نسبت به ایستگاه 
زاهدان کمتر است. علت اين موضوع را می‌توان در سری کامل و بدون انقطاع شاخص FDSD‏ در ایستگاه ایرانشهر 
نسبت به ایستگاه زاهدان جست‌وجو کرد؛ به طوری که در ایستگاه ایرانشهر تنها در دو مقطع زمانی تابستان و پاییز 
سال‌های ۲۰۰۵ و ۲۰۱۱ دو فصل متوالی بدون روزهای همراه با طوفان گردوغبار رخ داده است. 

در جدول Y‏ اعداد داخل پرانتز روبه‌روی polio‏ معیارهای ارزیابی رتبه‌بندی روش‌های مختلف پیش‌بینی فراوانی روزهای 
همراه با طوفان‌های گردوغبار را نشان می‌دهد. بر این مبناء در همة ایستگاه‌های مورد مطالعه به‌منظور پیش‌بینی شاخص 
۴02 روش سیستم استنباط عصبى-فازى تطبیقی (ANFIS)‏ بیشترین تعداد رتبة اول (بهترین حالت پیش‌بینی) را داشته 
اسان انهه بها :رز راتت Uu JUL DNE EC EE E CUTS‏ 
وجود خواهد داشت. بدین ترتیب مدل‌های سری زمانی dHolt-Winters SARIMA‏ و روش ۸8۴ اصلاً رتبة نخستی به 
خود اختصاص نداده‌اند؛ در حالی که روش ANFIS‏ ۲۰ بار در رتبة نخست قرار گرفت. درمحموع و با اين حساب» روش 
5 با اختلاف زیادی بهتر از روش‌های دیگر پیش‌بینی بوده و بالاترین دقت پیش‌بینی را داشته است. نکتة دیگری که 
بايد به أن اشاره کرد عملکرد مناسب شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ است. هرچند که اين روش در دو 
ایستگاه خاش و سراوان به‌ترتیب عملکرد ضعیف‌تر و lg‏ یکسانی با مدل سری زمانی Holt-Winters‏ داشته در غیاب 
روش 4۸۴15 با W‏ بار قرارگیری در رتبة دوم» بهترین عملکرد را به خود اختصاص داد. همچنین, نتایج دیگری از جدول ۴ 
قابل برداشت است که انتخاب مدل‌های ييجيدهتر به‌عنوان مدل بهينة پیش‌بینی کننده در ایستگاه‌های مورد مطالعه است؛ به 
نحوی که برای بيش بينى شاخص FDSD‏ در هم ایستگاه‌ها از مدل ۳ و L) Y‏ سه و چهار گام تأخیر) استفاده شده است. با 
توجه به اینکه بادهای ۱۲۰روزه عامل اصلی وقوع طوفان‌های گردوغبار در شرق و جنوب شرق کشورند» افزایش سرعت و 
E‏ یهام مزر رن tona‏ و ad dioe sus‏ درا ویک هرد 
قدرتمند برای وقوع طوفان‌های بعدی بوده كه در مطالعات فرج‌زاده اصل و علیزاده (۱۳۹۰)» راشکی و همکاران 4Y* W)‏ 
سبحانی و همکاران (Y)‏ یارمرادی و همکاران (۱۳۹۷)» عراقی‌نژاد و همکاران (WAY)‏ و قربانی و مدرس (WA)‏ بدان 
اشاره شده است. اين ذرات به‌جامانده از طوفان‌های قبلی می‌تواند دليل مهمی برای تأثير طوفان‌های گردوغباری چند فصل 
گذشته برای شکل‌گیری طوفان‌های گردوغبار در فصل آتی باشد. 

شکل ۸ polio‏ مشاهداتی و پیش‌بینی شاخص FDSD‏ را برای داده‌های آزمایش نشان می‌دهد. اين شکل بپبود 
عملکرد مدل‌های مورد استفاده به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار را با به کارگیری روش 
5 در همه ایستگاه‌ها به‌خوبی بیان می کند. به طوری که در همة ایستگاه‌های سینویتیک مورد مطالعه در استان 


سیستان و بلوچستان» روش ANFIS‏ بهترین مطابقت بين مقادیر پیش‌بینی و مشاهداتی شاخص FDSD‏ را داراست. 
به ۴۵ درجه (مدلهاى سرى زمانى SARIMA‏ و Holt-Winters‏ و روش (RBF‏ به یک خط دقيقاً ۵ درجه (روش 
(پیشروی از ایستگاه سراوان do‏ سمت ايستكاه زابل)» ارتباط بين مقادیر مشاهداتی 9 پیش‌بینی‌شده در همه روشها 


(مدل‌های سری زمانی و روش‌های هوش مصنوعی) مطابقت بیشتری با نيمساز ربع اول بيدا می‌کند. 
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مقايسة عملکرد مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ با روش‌های هوس مصنوعی در پیش‌بینی AAY e‏ 


مقايسة میانگین polio‏ مشاهداتی و پیش‌بینی‌شده بر اساس آزمون 2 (جدول ۵) نشان می‌دهد که فرض صفر مبنی 
بر برابربودن میانگین سری polio‏ مشاهداتی و ييش بينىشدة شاخص FDSD‏ در هیچ کدام از ایستگاه‌های مورد مطالعه 
بر اساس روش‌های ANFIS‏ و RBF‏ در سطح خطای Y‏ درصد و بر اساس مدل‌های سری زمانی SARIMA‏ و Holt-‏ 
۹ در سطح خطای ۵ درصد رد نمی‌شود. بنابراین» می‌توان نتیحه كرفت که don‏ روش‌های استفاده‌شده (به‌ویژه 
روش‌های ANFIS‏ و (RBF‏ به‌منظور پیش‌بینی شاخص FDSD‏ میانگین سری زمانى مشاهداتی را حفظ کرده‌اند (جدول 


۵ اين ویژگی می‌تواند تاکیدی بر کارایی روش‌های هوش مصنوعی در پیش‌بینی فراوانی روزهای همراه با طوفان 


گردوغبار باشد. 
جدول ۵. آزمون مقايسة میانکین سری زمانی مشاهداتی و پیش‌بینی‌شده شاخص FDSD‏ 
EUM ۲‏ آمارخ 2 
ایستگاه ‏ , ۳ Holt- Holt-‏ 
مشاهدا 5 ANFIS RBF SARIMA ANFIS 1 : SARIMA‏ 
Winters Winters‏ 
aha‏ د ۳۹۷ PA. PEE A‏ ۲ ۴ ۶ ۵ ۷ 
خاش ۴ EE‏ ان EE WEE ۳۹۷ mc iur‏ 
ایرانشهر YY ۰/۳۱۵ ۰۴١۱ EA «(NAA YYY ۰۴۵۸۵۵ EA UM‏ 
زاهدان ۵۰۲+ ARAN UY) YYY A ۱ ۰۵١ UNE ٠ «(NAA‏ 
زابل «NY YAY «YYY ۰۶۹۸ +۵۲ «AM‏ ۰*۵۱ ۰/۶۱۵ ۰۷۰۵ 
نتيجه کیری 


هدف از اين پژوهش مقايسة عملکرد مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ با روش‌های هوش مصنوعی 
شامل شبكة عصبی مبتنی بر توابع ib‏ شعاعی (RBF)‏ و سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ به‌منظور 
پیش‌بینی فراوانی روزهای همراه با طوفان‌های گردوغبار (FDSD)‏ در مقیاس فصلی در ينج ایستگاه سینوپتیک استان 
ks n‏ بكو زع و اه usu Ve EE‏ ایس فا al‏ کش یه اا EA‏ مود 
مطالعه روش ANFIS‏ بهترین عملکرد را نسبت به سایر روش‌ها داشت. همچنین. با افزايش شاخص FDSD‏ در 
ایستگاه‌های مورد مطالعه» دقت پیش‌بینی همه روش‌ها افزايش یافت. با اين توضیح که همه روش‌ها در ایستگاه زابل» که 
با ۷۱۱ روز در jb‏ زمانی ۲۵ساله (۲۰۱۴-۱۹۹۰) بیشترین فراوانی را به همراه داشت» بهترین عملکرد را از خود نشان 
دادند. غير از ایستگاه‌های خاش و سراوان» كه کمترین فراوانی روزهای همراه با طوفان گردوغبار را به‌ترتیب با مقادير 
۰ و ۱۰۶ روز در بازةٌ زمانی ۲۵ساله به خود اختصاص دادند. در دیگر ایستگاه‌های مورد مطالعه. روش‌های ۳15( 
Holt-Winters RBF‏ و SARIMA‏ به‌ترتیب رتبۂ اول تا چهارم را در پیش‌بینی شاخص FDSD‏ به‌دست آوردند. با کاهش 
فراوانی روزهای همراه با طوفان گردوغبار در ایستگاه‌های سراوان و خاش مدل سری زمانی cJ 94  Holt-Winters‏ 
عملکرد تقريباً یکسان و بالاتری نسبت به روش RBF‏ از خود نشان داد که حاکی از قابلیت بالای اين مدل به‌منظور 
پیش‌بینی مقادیر كم شاخص FDSD‏ بود. همچنین» نتایج بيانكر اين موضوع بود كه مدل سری زمانی SARIMA‏ نسبت 
به ساير روش‌های بيش بينى قابلیت بالایی در پیش‌بینی شاخص FDSD‏ در هیچ یک از ایستگاه‌های مورد مطالعه نداشت. 
همچنین, با وجود وقوع کمتر فراوانی روزهای همراه با طوفان گردوغبار در ایستگاه ایرانشهر نسبت به ایستگاه زاهدان» 
همه روش‌های پیش‌بینی شاخص FDSD‏ بر اساس معیارهای ارزیابی» عملکرد بهتر و دقیق‌تری نسبت به ایستگاه زاهدان 
داشتند که علت أن را می‌توان در وجود سری کامل و بدون انقطاع شاخص FDSD‏ در ایستگاه ایرانشهر جست‌وجو کرد؛ 


به طوری که در ایستگاه ایرانشهر فقط در دو مقطع زمانی تابستان و jek‏ سال‌های ۲۰۰۵ و ۲۰۱۱ دو فصل متوالی عاری 
از روزهای همراه با طوفان گردوغبار رخ داد. از طرفی» در پیش‌بینی شاخص FDSD‏ در استان سیستان و بلوچستان» مدل 
بهينة پیش‌بینی کننده مدل پیچیده‌ای بوده است. به صورتی که برای همة ایستگاه‌های مورد مطالعه. مدلی که از سه يا 
کار گام eb‏ هر ن ےم استفاده کرد هعتوان رین مدل فش كدو od os‏ از این رف دراك به‌عامانده از 
طوفان‌های قبلی می‌تواند دلیل مهمی بر تأثير طوفان‌های چند فصل گذشته بر شکل‌گیری طوفان‌های گردوغبار در 
فصل‌های آینده باشد. يعنى طوفان شکل گرفته‌شده در حال حاضر ممکن است به خاطر ذرات رسوب کرده و ته‌نشین‌شده 
از سه يا چهار فصل گذشته باشد که اکنون به خیزش درآمده است. همین دلیل علمی که اکنون بیان شد در دیگر 
تاش ان انه dala‏ میا hacc Sulis dafs‏ اه Bebidas. 5555 Scl‏ 
سنجش صحیح عملکرد مدل‌ها در پیش‌بینی شده هم علت شکل‌گیری طوفان گردوغبار در مناطق مستعد اين پدیده را 
آشکار کرد؛ زيرا ذرات به‌جامانده از طوفان‌های گردوغبار قبلی خود عاملى بر شکل‌گیری طوفان‌های گردوغبار فعلی و آتی 


.+ 
هسنند. 
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